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INTRODUCTION 

This report presents a method for generating all the orbital parameters 
needed for the thermal analysis of spacecraft in Earth monitoring 
Sun -synchronous orbits. The level of complexity is suitable for initial 
planning of missions. The mechanics of Sun- synchronous orbits and their 
impact on the thermal design of space experiments is discussed in detail. A 
Sun- synchronous Qrbit Analysis £rogram (SOAP) was developed to calculate the 
orbital parameters necessary for a complete thermal analysis. An example case 
run with this program is discussed. 

ANALYSIS OF SUN-SYNCHRONOUS ORBITS 


An overview of Sun- synchronous orbits is presented here with particular 
emphasis on their application to Earth-monitoring missions. Additional 
analysis of Sun- synchronous orbits can be found in reference 1. A 
Sun- synchronous orbit is one whose right ascension of ascending node precesses 
at a rate which matches the average apparent rotation rate of the Sun about 
the Earth. This type of orbit is desirable for many types of Earth-monitoring 
missions because it allows for near-global coverage of the Earth's surface and 
atmosphere and repetitive latitude -longitude observations at near- constant 
conditions of solar illumination. Also of special significance for Earth 
observing satellites, are the times during which the Sun is rising or setting 
on the Earth's horizon, as viewed from the spacecraft. During sunset and 
sunrise it is possible to obtain vertical profiles of many different 
atmospheric constituents which attenuate the Sun's radiation. A spacecraft in 
Sun- synchronous orbit, with a typical low Beta angle 1 (angle between the solar 
vector and the orbit plane) , will continue to experience a sunrise and sunset 
every revolution, for as long as the orbit is maintained. In the case of a 
purely Earth monitoring experiment, atmospheric data can be recorded over the 
same area under the same lighting conditions at a constant time increment. 


Sun- synchronous orbits are retrograde, that is, the velocity component of the 
spacecraft in the equatorial plane opposes the rotational motion of the Earth. 
Thus, the orbit inclination (angle between the orbit plane and equatorial 
plane) will always be greater than 90 degrees. For circular orbits, the 
inclination is only a function of the orbit altitude. The variation of 
Sun- synchronous orbit inclination for circular orbits up to an altitude of 
1600 kilometers is shown in Figure 1. 


■^The Beta angle defined in this paper will always refer to the angLe between 
the orbit plane and the Sun, Other publications often use the complement to 
this angle. 
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Another special characteristic of a Sun-s /nchronous orbit is that the orbit 
repetition factor (Q) is identical to the number of orbit revolutions in one 
day. Therefore, assuming a circular orbit, any Sun- synchronous orbit has the 

special property that its ground track wiLl repeat precisely at intervals of a 

whole number of days. Table 1 shows the number of orbital revolutions per day 
as a function of orbit inclination or orbit altitude. This quantity Q 

is important for defining the longitude and latitude ground coverage of an 

Earth monitoring satellite. For example, an integer value for Q means that a 
particular ground track will be duplicated the next day. If Q includes a 
fraction, such as one-fifth, then the initial longitude coverage will be 
duplicated every five days. Therefore, long repeat cycles will allow for 
good longitude coverage, but poor repeat coverage. Similarly, short repeat 
cycles give good repetition, but poor longitude coverage (ref. 2). 

It is necessary to note that, though a Sun -synchronous Earth monitoring orbit 
implies a fixed geometry with respect to the Sun, the solar radiation flux 
input is not necessarily fixed, but actually varies considerably. This is due 
to the satellite's path which crosses around the Earth and maintains a 
constant view of the Earth on one side. Though no single side of the 
satellite is fixed with respect to the Sun, there is often one sida that 
receives a high average solar flux during an orbit. For example, if the 
satellite orbit plane is in line with the Sun, the satellite will receive 
direct solar radiation on its zenith side for the part of the orbit when the 
satellite is between the Earth and the Sun. As the satellite moves past the 
polar region it will go behind the Earth, thus being blocked from the Sun's 
radiation entirely. This type of orbit configuration can cause problems for 
optical experiments using nadir detectors with a horizon to horizon view. 

In this case, the nadir, or earth facing side of the experiment, will 
experience direct solar flux during sunrise and sunset as seen from the 
satellite as it passes around the Earth. The other extreme case occurs when 
the satellite orbit plane is perpendicular to the Sun. This will cause one 
side of the satellite to receive full solar flux. Similarly, the opposite 
side of the satellite would be facing deep space during the entire orbit. 

As with any space experiment, the solar illumination or radiation flux input 
has a considerable Impact on the thermal and mechanical design. Therefore, it 
Is appropriate to discuss the orbital alignment with the Sun, specifically as 
it applies to an Earth monitoring Sun- synchronous experiment. In chis case, 
the initial timing of the orbit is important, since it will determine the 
envelope of angles between the solar vector and the orbit plane. The initial 
timing of a Sun- synchronous orbit is often referred to as the "orbit time", 
which represents the longitudinal difference between the sun and the point 
where the spacecraft crosses the equator (15 degrees per hour). For example, 
if an orbit time of local clock noon is chosen, then the orbit plane would be 
closely In line with the Sun and the angle between the Sun and the orbit plane 
would be very small, even approaching zero at times. Each time the satellite 
passes the equator the local clock time at this crossing position on the Earth 
would be noon. Similarly, an orbit time six hours later would produce angles 
near 90 degrees and an envelope that generally Includes larger angles. This 
type of orbit could create thermal problems since a particular side of the 
experiment would be facing the Sun at all times, providing the nadir side 
maintains a constant Earth view. In any event, the characteristics of a 
Sun- synchronous orbit have a direct effect on the design of a space 
experiment . 
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It Is customary that the orbit time be expressed as either an ascending node 
or descending node. This describes the path of the satellite as it crosses 
the equator. In the case of an ascending node orbit, the satellite passes the 
equator traveling from South to North. Similarly, a descending node orbit 
crosses the equator traveling from North to South, Thus a 2:00 pm ascending 
node orbit is equivalent to a 2 : 00 am descending node orbit. 

SUN-SYNCHRONOUS ORBIT ANALYSIS PROGRAM 

A Sun- synchronous Orbit Analysis Program (SOAP), written in FORTRAN language, 
was created to calculate orbital parameters for an entire year for a fixed 
orbit altitude and orbit time. The source program has been written in the 
form of a subroutine, designed to be called from any other computer code. 

The subroutine and a sample calling program are supplied in Appendix A and an 
example output is presented in Appendix B. The input and output calling 
arguments, for the SOAP subroutine, are shown in Appendix C and the governing 
equations used to calculate orbital parameters are presented in Appendix D. 

The program is valid for Sun -synchronous orbits with launch dates after 1950 
including leap years. Inputs for the program include the year of the launch, 
experiment altitude and Sun- synchronous orbit time. Output data includes 
orbit inclination angle, time in hours for a complete orbit or revolution 
time, a set of data for the maximum and minimum Beta angle, and orbital data 
for the entire year including solar declination, right ascension of the Sun, 
and longitude of the ascending node. The orbital parameters calculated by 
this program are extremely tedious to do by hand and require repetitive 
calculations for every day of the year. Thus, this program provides an 
efficient method for determining all of the orbital parameters needed for a 
thermal radiation analysis. 

It was stated previously that a Sun- synchronous orbit allowed for a fixed 
inclination angle and a fixed revolution time. All of the other parameters 
are a function of orbital position and the day of the year. A typical set of 
output data will contain a maximum and minimum Beta angle for the year of 
interest. Though the year and altitude have a small effect on the Beta 
angles, the Sun- synchronous orbit time has the largest effect on t lese angles 
and their corresponding magnitudes. For example, a 6:00 orbit produces a 
large variation in the Beta angle with the average angles being high, even 
approaching 90 degrees in some cases. An orbit time of 12:00 would produce 
the smallest variation in Beta angle and the maximum and minimum angles would 
both be below 10 degrees. The yearly variation in Beta angle Is due to 
seasonal variation of the Earth's position with respect to the Sun. In the 
case where the orbit plane is inclined at a high angle with respect to the Sun 
(orbit times near noon), the Earth's seasonal tilt will create large 
variations in Beta angle. Similarly, when the orbit plane is closely aligned 
with the Sun (orbit times near 6:00), the Earth's tilt has little effect on 
the Beta angle. Figure 2 shows the variation in Beta angle for different 
Sun- synchronous orbit times. These results, calculated by the SOAP computer 
program, are for a fixed altitude and year. 


^The source code can be acquired on either diskette or magnetic tape by 
contacting the author of this publication. Note, the author will not be 
responsible for any changes in the source code made by the user. 
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THERMAL ANALYSIS APPLICATION 


A specific example case was chosen to illustrate the use of the 
Sun- synchronous Orbit Analysis Program (SOAP) in calculating the orbital 
fluxes incident on a Sun -synchronous space experiment. The example utilizes a 
radiation analysis program to predict heat: fluxes on a cube for one complete 
orbit. Results from the orbital analysis program, described in the previous 
section, are used as input to this radiation program. An example orbit will 
be presented and the orbital heating data will be summarized. 

The orbit chosen for this application was a 1:30 p.m. ascending node 
Sun- synchronous circular orbit with an altitude of 705 km and a launch date in 
1998. Representative values for the solar constant, Earth infrared radiation 
and the Earth albedo were used. These values are commonly used by thermal 
analysts for spacecraft design. It was also assumed that the spacecraft 
surface was black, such that the emlssivity and absorptivity were 1.0. This 
orbit data was used in the SOAP computer code to calculate pertinent orbit 
parameters required by the Thermal Radiation Analyzer System (TRASYS) program 
(ref. 3). These parameters are summarized in Table 2. The TRASYS program was 
employed to define the incident infrared, solar and Earth albedo heat fluxes 
acting on a space experiment for one complete orbit. Thermal flux data were 
generated for one orbit on two distinct days of the year when the Sun's 
location, relative to the experiment, yields the highest and lowest Beta 
angles. These two days represent the highest and lowest average heat fluxes 
for the year. The ultraviolet (UV) , which includes the solar and albedo 
fluxes, and the Earth infrared (IR) flux contributions are shown in Table 3. 
These heat fluxes can be used as the basis for further thermal analysis 
including sizing radiators or predicting instrument operating temperatures. 
Since these two sets of data contain the highest and lowest average heat 
fluxes for the year, the thermal analyst can be assured that all orbital 
heating fluxes will fall within these limits. 

CONCLUSIONS 

The fundamentals of an Earth monitoring Sun- synchronous orbit were 
presented. A Sun- synchronous Orbit Analysis Program (SOAP) was developed to 
calculate orbital parameters for an entire year. The output from this program 
provides the required input data for the TRASYS thermal radiation computer 
code, which in turn computes the Infrared, solar and Earth albedo heat fluxes 
incident on a space experiment. Direct incident heat fluxes can be used as 
input to a generalized thermal analyzer program to size radiators and predict 
Instrument operating temperatures. The SOAP computer code and the thermal 
analysis methodology presented in this paper should prove useful to the 
thermal engineer during the design phases of future Earth monitoring 
Sun -synchronous space experiments. 
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REPETITION FACTOR (Q) 

ALTITUDE (KM) 

INCLINATION ANGI E (DEG) 

12 

1684.6 

102.94 

121/8 

1629.1 

102.63 

121/4 

1574.4 

102.33 

121/2 

1467.9 

101.75 

13 

1265.3 

100.71 

13 1/8 

1216.6 

100.47 

13 1/4 

1168.7 

100.24 

13 1/2 

1075.2 

99.79 

14 

896.5 

98.99 

14 1/8 

853.4 

98.80 

14 1/4 

811.0 

98.62 

14 1/2 

728.1 

98.28 

15 

569.1 

97.64 

151/8 

530.8 

97.49 

151/4 

492.9 

97.35 

15 1/2 

418.8 

97.07 

16 

276.3 

96.57 


Table 1. Circular Sun-Synchronous Orbit Parameters for Specific Repetition Factors 


PARAMETER 

DATE 


July 9, 1998 

November 5, 1 998 

Days After January 1 

190 

309 

Longitude of Ascending Node 

130.77 

248.07 

Right Ascension of the Sun 

108.19 

220.10 

Solar Declination 

22.39 

-15.60 


Orbital Constants: 

Orbit Altitude = 705 km 

Sun-synchronous Orbit Time = 1 :30 pm Ascending Node 
Revolution Time = 1 .65 hours 
Eccentricity = 0.0 

Thermal Analysis Constants: 

Solar Constant = 1352 W/m*m 
Earth Infrared Radiation = 236.6 W/m‘m 
Earth Albedo = 0.30 

Incident Flux Properties: 

Absorptivity = 1 .0 
Emissivity = 1 .0 

Table 2. TRASYS Input Parameters 
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Table 3. Experiment Heat Flux Data 
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Appendix A 

SOAP Computer Program Listing 


Sun- Synchronous orbit Analysis broom (SO*#) 
written In the torn ot a FORTRAN SUBROUTINE 

Written by: Or lan KH lough - NASA Langley Research Center 

Mil Stop ot 

Heap ton. va (464) **4-7047 

This program calculates all of the necessary orbital parameters 
lor a sun- synchronous space experiment, input for the prograa 
includes the year, experiment altitude, and sun- synchronous orbit 
time. Output for the prograa includes the BETA angle (Angle 
between the sun and the experiment orbit plane) for each day of 
the year, and the Mxlaum and minimum BETA Angles for that year. 
These maximum and minimum SETA Angles are used to calculate 
orbital parameters such as: Julian Date. Days Relative to January 
i, solar Declination. Orbit inclination, longitude of the Ascending 
Node, tight Ascension of the Sun. and the Orbit Time. These 
parameters can be used to run any radiation analysis program that 
calculates radiation fluxes over one orbit, on a given day. 


PROGRAM MIN 

• THIS IS AN EXAMPLE CALLING PROGRAM 


WRITE(2.2BO) 00X73(2) 

2to format (IX.' Days Relative to January i ■ u) 

WRtTE(2.2*0> OfCL I (COUNTS! 2) ) 

2*0 format ux, * Solar Oocllnatlon.deg. « \F4.2) 

WRITEf 2.360) LONCAS (C0UNTS( 2) ) 

300 FORMT (IX,' Longitude of Ascending Node. deg, * ’.F0.2) 

WRITE (3, 310) RASSLMK COUNTS! 2) ) 

310 FORMAT (ix. ' Right Ascension of the Sun, deg. ■ .Fa. 2 ) 

WRITEC2, 320) 

320 FORMT { IX. ’* ••••••• •••••••• * 

WRITE! 2 . 340) 

340 FORMT (lx. 'Day of Year '.2X. Beta Anglo', 2X. 'Solar Declination’ 
>2X. 'Right Ascension’ ,2X. ' Longitude of*./.' *.2x, 

,2X. ,4X. 'of the Sun ’.ax. 

>' Asctndlng Node’ , /.4Sx, . ix. ’ «■»■•••■•■«■■') 

00 400 UI.MSwLEAP 

WRITE(2,350) I . SETa( i ) , DCCl 1(1). RASSUN( I ) . LONGA Si I ) 

3 SO FORMT (3X.I4.tX.F*.2, 10X.FB.2, I2X.FR.2, I1X.FR.2) 

400 CONTIMJE 

CLOSE (unit. 2) 

STOP 

EAO 

SUBROUTINE SOAPiVlAR. ALT. TLAINCH. LEAP. MONTHS. DAYS. COUNTS. I NCI. 

> OTlMC.BETA.OfO. I, LONCAS. RASSUN. (ERROR) 


DIMENSION 6CTA<364}.DECLK34*}.10NGAS(34*),RASSUN(344) 
DIMENSION MONTHS! 2). DA VS( 2) .COUNTS! 2) 


REAL iETA.DECll . LONCAS . RASSUN 
REAL ALT.TLALNCH, INCL .OTIME 


INTEGER LEAP. MONTHS. OAYS.COtNTS. YEAR 


CHARACTER *20 OUTFILC 


I CONTINUE 


PRINT 10 

10 FORmT (/.' sun- Synchronous Orbit Analysis Program (SOaP)'./. 

> ' 

PRINT 20 

20 format* ix. enter an output file name In sincie quotes:’) 

READ ‘.OUTFILE 

OPEN! INI TO . F I L E «OUTF i L E ) 

PRINT 30 

30 FORmT! ix. 'Enter the Year:') 

READ • .YEAR 
PRINT 40 

40 FORMT! IX, 'Enter the Experiment Altitude 
READ * . ALT 
PRINT SO 

SO FORMT! ix. ‘Enter the Sun- Synchronous orbit Time:’) 

PRINT *0 

to FORMT! IX. Example » 13.5 {* 130 p« a sc ending Node)’) 

RCAO vtlaunch 

CALL SOAPtVEAR. ALT. TlAUNCM.lCAP. MONTHS. DAYS. COUNTS. INCL , OTIME . 
> BETA.DECL I . LONCAS. RASSLN, I ERROR) 

IF (I ERROR NE. 0 ) GO TO 1 


WRITE (2.70) 

70 FORMT (IX, ’ 

WRlT(<2,aO) 

SO FORmT (ix. ' Sun- Synchronous orbit Analysis Program (SOAP)’) 
WRIT|(2.aO) 

90 FORMT * 


WRITE! 2, 100) YEAR 

100 FORMT (IX, ' Year of Launch ■ .14) 

WRITE f 2. 120) ALT 

120 FORMT (IX.* Altitude (Rm) « '.FS.1) 

WRITE(2. 140) HAUNCH 

mo FORmT (ix, ‘ Sun-Synchronous orbit rime « *.FS.2) 

WRITE(2. ISO) 

ISO FORMT (IX. '•••••••••■■•• 

Wt!TE(2.!tO) INCL 

t*0 FORMT (ix. * Orbit incl Inat Ion Angle. deg. * \F5.2) 

WR1Tf(2, 170) OTIME 

170 FORMT (IX.' Orbit Time (hours) ■ \F4.2> 

WRITE (2, 1*0) 

1*0 FORMT (IX. 

WRITE (2. 190) *ETA(COUNTS{ I ) ) 

190 FORmT { »X. ' Maximum Beta Angle. deg. * \F4.2> 

WR 1 Tf ( 2 200) MONTHS!! > DAYS! I ) 

200 FORmT (tx.' Date of Maximum Beta Angie • Ml. /Ml) 

WRITE! 2. 210} COUNTS! 1 ) 

210 FORMT (ix. ' Days Relative to January I • \I4> 

WR(TE(3. 220) OCCL I (COUNTS! t ) } 

220 FORMT (ix.* Solar Dec I inat ion. deg. • ’ . fd2> 

WRITE (2. 230) LONCAS (COUNTS! I}) 

230 FORmT (IX,* Longitude of Ascending Node. deg. ■ .F0.2) 

WRITE (2, 240) R ASSUfl COUNTS! U) 

240 formT (ix.' Right Ascension of the Sun. deg. ■ \F*.2) 

WRITE (2. 290) 

250 FORMT (IX. ) 

WRITE (2 . 200) MTA( COUNTS! 2) ) 

200 FORMT (IX,* Rlnlmum Beta Anole.Oeg. ■ .FB I) 

WRITE! 2 270) M94THS!2) DAYS! 2) 

270 FORMT (ix. * Oett of Minims* Rota Angle ■ Ml. / ’. 12) 


Sun-Synchronous Orbit Analysis Program (SOAP) 
written In the form of a fortran surroutinc 

written by: Brian xl Hough - NASA Langley Research center 

Mail Stop 431 

Hampton. VA 33445 (B04) *44-7047 

This program calculates all of the necessary orbital parameters 
for a sun- synchronous space experiment, input for the program 
Includes the year, experiment altitude, and sun- synchronous orbit 
time. Output for the program includes the BETA Angie (Angle 
between the sun and the experiment orbit plane) for each day of 
the year, and the maximum and minimum beta Angles for that year. 
These maximum and minimum BETA Angles are used to calculate 
orbital parameters such as: Juiianoate. Days Relative to jamary 
l. Solar Declination, orbit inclination. Longitude of the Ascending 
Node. Right Ascension of the Sun. and the Orbit Time These 
parameters can be used to run any radiation analysis program that 
calculates radiation fluxes over one orbit, on a given day. 


• VARIABLE DEFINITION: 


• INPUT: 

YEAR « YEAR REQUEST >• 1950 

alt * Experiment Altltudg 

tlalnch « Sun-Synchronous orbit Time (Hourj past Midnight) 

• aumn: 

LEAP -I IF LEAP VEAR; 0. IF NOT LEAP YEAR 
MONTHS - MONTH INDICATOR FOR MAX. MIN BETA 
DAYS v DAY INDICATOR FOR MX. MIN BETA 
COUNTS * INDEX INDICATOR FOR MAX. MIN IN ARRAYS 
INCL « inclination Angle 

OTImE • Orbit time (hours) 

BETA v BETA Angie (Angle between the Sun and the Orbit Plane) 
DCCli ■ Solar Declination 
loncas ■ longitude of the Ascending Node 
rassln ■ Right Ascension of the sun 

I ERROR • 0. IF NO ERROR: 1. ERROR IF DATE < 1950 CAL. NOT VALID 


VARIABLES USED IN LOGIC: 
jldate - Julian Date 
mlonG ■ Man Longitude of the Sun 

manor ■ mean Anomaly 

ECLONC ■ fcl iptic Longitude 
041 EC ■ obliquity ol Ecliptic 

racreEn » Right Ascension of the Greenwich meridian at Midnight 
xnorm • vector Normal to the Experiment (Ra-oecl coord. System) 
sln ■ Solar vector (ra-occl coord, system) 


DIMENSION BETA! 344) DEO. I ( 344) . LONCaS! 344) . RASSUN! 344) 
DIMENSION COUNTS! 2) .MONTHS! 2) , OAYSf 2) 


* MX I MUR THEN MINIMUM STORED IN ARRAYS 

INTEGER VE AR. LEAP, MONTHS. DAYS, C0LNT3. (ERROR 


REAL ALT. TIAUNCH. INCL, OTIME 
REAL BETA.DECL I .LONCAS, RASSUN 


REAL Ml ONC. MANOR, LESS. INC1. INC2. 

» P I. RAO. ARC, ECLONC. JUDAT. JLDATE. NUNMY. 

> ORLfC.lASI . RAS2.RAGREEN.ANC1 . ANC2 . XN0KM I , XN0RM2 . XNQRM3 . 

» XSUN. VSLN.ZSUN. XN9UN, BETAMAX. BCTAMIN 

INTEGER COUNT. COUNT 1. C0UNT2. 1 . NBIRST.NL AST 
INTEGER DAVOtntJ) 
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DATA DAVCNT/0.31.39.90.120. I S f . Ill , 213,243. 273. 304 , 334 . 366/ 


• 0CF1NR JU.OAT • 12:00 NOON JAN. 1. 1050 

m 

DATA JULDAT/243J282 3/ 


PRINT 2 

2 FORMAT</. IX. '•**•••••• PLEASE WAIT - PROCESSING ./) 

PI *3 . 14 150265339 
LESS-O. 

LEAP-0 
RAO* 100. /PI 


• CHECK YEAR REQUEST • 


I ERROR * 0 

IHYEAR.lt. 1950) THEN 

PRINT YEAR REQUEST LESS THAN 1950 . . CANNOT \ 
> 'CALCULATE' 

I ERROR * I 
GO TO 200 
ENDIF 


• CALCULATE INCLINATION ANCLE AND TIME 


INC 1 * RAD*ACOS(0 . 9050473* (6378 . 143+ALT) * *3.5 

> / ( -2004 .3874E+11)) 

INC2- 1 ♦< 66058 . 33/( (ALT+6378 . 145) • • 2) J • { 1- 1 . 3*<SIN< INC1 

> /RAO)* *2)) 

INCL * RAO* ACOS( (COS( INC 1 /RAO) )/ InC2) 

OT I ME -0.0 

IF {(INCL LT. 97. 641). AM). { INCL. GE. 96. 566)) OTlMf* 

> 24/< < INCL -97 . 64 1 ) / ( 94 . 566-97 64 1 ) ♦ 1 3 ) 

IF { ( INCL .LT. 96 . 966) .ANO, ( INCL .GE . 97 . 641 ) ) OTlMf* 

> 24/( ( INCL -96. 988)/(97. 64 1-98.988)+ 14) 

IF ( ( INCL . LT. 100.706) .ANO. ( INCL .GE . 96 . 966) ) OTIME- 

> 24/U INCL- 100. 706) /(9B.9BB- 100. 706)* 13) 

IF (( INCL. LE. 102.944) AND. ( INCL CE . 100.706)) OTIME- 

> 24/ ( ( INCL -102 . 944)/( 100.706-102. 944)- 12) 


IF ( (ANC2 . Gf . 90. 0) .AM). (ANC3.LT. 140.0)) ANG2-ANG7-90 . 

IF ((ANC2.CC. IU.0) .AND. (ANC2.LT. 270.0)) ANC2*ANC2- 1 80 . 
ANG2* ANC2- 270 . 

END IF 

IF ( (ANCI . GC .90.0) .AND. (ANGI . LT. 160.0)) THEN 

IF ( (ANC2 . Cf . 0. 0) .AM). (ANC2.LT.9O.0)) ANC2-ANC2+90 . 

IF ((ANG2. CE. 90.0) AND. (ANC2.LT. 160.)) COTO 140 
IF (<ANC2.CE. 180.0). AND. (ANC2.LT. 270.)) ANC2-ANC2-90 . 
ANG2-ANG2- 1 60 . 

ENDIF 

IF ( (ANCI .CE . 160. ) .AND. (ANCt LT. 270. )) THEN 

IF ((ANC2.GE .0.0) AND. (ANC2.LT. 90.0)) ANC2-ANC2+ 1 80 . 

IF ( (ANC2.CE. 90.0) AND. (ANC2.LT. 180.0)) ANC2-ANC2+90 . 

IF ( (ANC2.CE . 180.0) . AND. (anC2.LT. 270.0) ) COTO 140 
ANG2-ANC2-90 . 

ENDIF 

IF ((ANCI. CE. 270) AND. (ANCI. IT. 360.)) THEN 

IF ((ANC2.CE.0.0) AND. (ANC2.LT. 90.0)) ANC2-ANC2+270 . 

IF ((ANC2.GC.90 0) .ANO.UNC2.1T. 180.0)) ANC2-ANG2 + 1 80 . 
IF ((ANC2.CE. 180.0) AND. (ANC2.LT. 270.0)) ANC2-ANC2+90 . 
ENDIF 

140 RASSUNC COUNT)* ANC2 


• BETA ANCLE CALCULATION 


XNORM1 -5 IN( INCL /RAO) *S IN( LONGaS ( COUNT ) /R AO) 

XNORM2-- SIN( INCL/RAO) 1 COS (L0NCAS( COUNT) /RAO) 
WORM3*COS( INCL /RAO) 

XSUN-COS ( ECLONC/RAO) 

YSUNMS I N( ECLONC/R AO ) *COS < OOL EC/R AO) 

ZSL**S I N( ECLONC/RAO) * S I N{ Oftt EC/R AO) 

WOINbRAD* ACOS ( WORM 1 • XSUN+ XNORM2' * YSUN+ WORMS- * ZSLN ) 
BETA (COUNT)* AES ( 90-XNSUN) 

150 CONTINUE 


• MAXIMUM - MINIMUM BETA ANGLE5 * 


BETAMAX*BETA(1) 
BETAM IN«BETA( 1) 


■ CALCULATE JULIAN DATE FROM JAN 1 1930 


NLUfDAY • (YEAR - 1930) * 363 . ♦ JULDAT ♦ INT< FIOATWAR- 1950 ) /4 . 
> ♦0.49) 


IF <M00< (YEAR-1930). 4) .EQ. 2) LEAP*1 


■ ORBITAL PARAMETER CALCULATION * 


00 ISO COUNT- 1 . 363-LEAP 
JLDATE * COUNT ♦ NLMDAY 
ARC* -4018.5 ♦ COUNT 

10 AMONG* 280.460+0. 9856474 *ARG+L ESS* 360 . 

IF (MLONC.CTO.O) COTO 20 
LESS*LESS+1 . 

COTO 10 
20 Lf SS*0 . 

30 MANOM* 357.528+0 9136003* ARC+L ESS* 360. 

IF (MANOM. CT. 0.0) GOTO 40 
LESS-LESS+1 
GOTO 30 
40 LE5S*0 . 

ECLONC* MLONC+ 1 . 91 3*S IN(MANOM/RAD) +0 , 02* S IN( 2* MANOM /RAO) 
OELEC* 23.439-0. 0000004* ARC 

OE CL I ( COUNT ) * RAD* A5IN( S IN{ ECLONC/RAO) *SIN(OELEC/RaO) > 
RASI* ( JLDATE - 24 1 5020. 0)/3652 5* 36000. 7689 
RAS2* 0.00038708* ( { JLDATE -24 15020 . 0) Z36525) **2 
RACREEN* 99 . 6909833+RAS 1+RAS2-LESS* 360 . 

IF (RACREEN LT. 360 0} COTO 60 
LESS-LESS+1. 

COTO 50 
IESS*0. 

LONCASt COUNT)* R AG* EEN+TLALNCH* 60. *0.25068447 -LESS *360. 
IF (LONCAS( COUNT) LT. 360.0) GOTO 80 
LESS-LESS+l 
GOTO 70 
LESS-O. 

RAS$UN( COUNT)* RAO* ATAN( COS (OBL EC/RAO) • TAN( ECLONC/RAO) > 
IF ((RASSUN( COUNT) .CT.O.O) .AN>. (RASSLN( COUNT) LT. 360.0)) 
GOTO 130 

IF (RASSUN< COUNT) .CT.360 . 0) COTO 100 

IF (RA5SUSI( COUNT) LT.O.O) RASSUN(COUNT).RASSUN(COUNT) ♦ 

LESS* 360 

LES5*LESS+t. 

GOTO 90 
LESS- I . 

RASSUN(COUNT) *RA5SUN( COUNT) -L ESS* 360 
LiSS-LESSd. 

IF (RaS 3LN( COUNT) . GT. 360.0) GOTO 120 
ANCI -ECLONC 
ANC2«i ASSUME COUNT) 

IF ((ANCI Cf .0.0) AM). (ANCt LT 90 0)> T>*N 

IF ( (ANC2.CE .0.0) .AM>. ( ANC2.LT . 90. 0) ) COTO 140 


30 

60 

70 

80 

90 

100 

120 

130 


00 >60 1-2. 365+lCAP 

IF (BETAMlN.CT. BET a{ I ) ) THEN 
C0UNT2- I 
BETAMIN-BETA(I ) 

ENDIF 

IF (BETAMAX.LT BETA( I ) ) THEN 
COUNT 1- I 
BETAMAX-BETa{ I) 

1*0 IF 
160 CONTINUE 


' FIND DATE OF MAX AND MIN • 


NFIRST - DAYCNT ( 1 ) 

00 170 1*1.12 
NLA5T ■ DAYCNT ( I ♦ 1 ) 
t FC I CT. I) THEN 
NL AST - NLAST + LEAP 
ENDIF 

IF(C0UNT1.CT. NFIRST .AND. COUNT! LE. NLAST) THEN 
COUNTS(I) - COUNT I 
DAYS ( 1 ) - C0UNT1 - NFIRST 
MONTHS ( t ) • I 
ENDIF 

IFIC0UNT2 CT. NFIRST AND. C0UNT2 IE NLAST) THEN 
C0UNTS{2) ■ C0UNT2 
0AVS(2) - COUNT 2 - NFIRST 
MONTHS( 2) ■ I 
ENDIF 

NFIRST • NLAST 
170 CONTINUE 

200 CONTINUE 
RETURN 
END 
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Appendix B 
SOAP Example Output 

** Sun- Synchronous Orbit Analysis Program (SOAP) ** 
•irkickicfck1rk*irlcirkirk*1ckick*1ck1rk*1tick*r k 1cic kirk 1ciWrirkTkirkick1r/cirk 

Year of Launch - 1998 
Altitude (km) - 705.0 
Sun- Synchronous Orbit Time - 13.50 

Orbit Inclination Angle, deg. * 98.21 
Orbit Time (hours) - 1.65 

Maximum Beta Angle, deg. - 29.04 

Date of Maximum Beta Angle - 11/5 

Days Relative to January 1 - 309 

Solar Declination, deg. - -15.60 

Longitude of Ascending Node, deg. - 248.07 

Right Ascension of the Sun, deg. - 220.10 

Minimum Beta Angle, deg. - 17.28 

Date of Minimum Beta Angle — 7/9 

Days Relative to January 1 - 190 

Solar Declination, deg. - 22.39 

Longitude of Ascending Node, deg. - 130.77 

Right Ascension of the Sun, deg. - 108.19 

*^oV^^*****-*r*-)V****************VSo!n!r>k“*^Vr>t************ 



Longitude of 
Ascending Node 

304.48 

305.47 

306.46 

307.44 

308.43 

309.41 

310.40 

311.39 

312.37 

313.36 


358 

25.24 

-23.42 

272.34 

296.36 

359 

25.13 

-23.40 

273.45 

297.35 

360 

25.02 

-23.37 

274.56 

298.33 

361 

24.91 

-23.34 

275.67 

299.32 

362 

24.79 

-23.29 

276.78 

300.30 

363 

24.68 

-23.24 

277.89 

301.29 

364 

24.57 

-23.18 

279.00 

3)2.28 

365 

24.46 

-23.12 

280.10 

303.26 
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APPENDIX C 


Input and Output Calling Arguments for the SOAP Subroutine 

SOAP has been written in FORTRAN in the form of a subroutine, designed to be 
called from any other computer code. The input and output calling arguments 
are described below in detail. 

USE: CALL SOAP (YEAR , ALT , TLAUNCH , LEAP , MONTHS , DAYS , COUNTS , 

INCL , OTIME , BETA , DECLI , LONGAS , RASSUN , I ERROR) 

INPUT: YEAR An input integer, designating the year, must be >“ 1950. 

ALT An input real variable, experiment altitude, km. 

TLAUNCH An input real variable, Sun- synchronous Orbit time (hours 
past midnight). i.e. 13.5 ( - 1 : 30 pm ascending node) 

OUTPUT: LEAP Output integer variable, - 0, YEAR is not a leap year; - 1, 

YEAR is a leap year. 

MONTHS Output integer array, dimensioned 2, containing the month 
indicator for maximum followed by the minimum BETA 
occurrence . 

DAYS Output integer array, dimensioned 2, containing the day 

indicator of the value MONTHS for the maximum followed by 
the minimum BETA occurrence . 

COUNTS Output Integer array, dimensioned 2, containing the Index 
pointers in the BETA, DECLI, LONGAS, RASSUN arrays for the 
maximum followed by the minimum BETA occurrence. 

INCL Output real variable, the inclination angle, deg. 

OTIME Output real variable, Orbit time, hours. 

BETA Output real array, dimensioned 366, angle between the Sun 

and the Orbit Plane, deg. 

DECLI Output real array, dimensioned 366, solar declination 
angle, deg. 

LONGAS Output real array, dimensioned 366, longitude of the 
ascending node, deg. 

RASSUN Output real array, dimensioned 366, right Ascension of the 
Sun, deg. 

IERROR Output integer, * 0, if no error; - 1, error if YEAR < 

1950. 
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APPENDIX D 


Governing Equations for the Calculation of Orbital Parameters 

Assuming a fictitious spherical Earth and a nearly circular orbit, the 
altitude of the spacecraft can be expressed as 

a - h + 6378.145 (Dl) 

where a is the altitude above the center of the Earth and h is the nominal 
altitude of the satellite above the Earth's surface in kilometers. Using an 
iterative approach, the orbit inclination can be found from the following 
equations (ref. 4) 

cl - cos' 1 ( 0.9856473 a 3 - 5 / -2.0645874 x 10 14 ) (D2) 

c2 - 1 + 66058.33 (1-1.5 sin 2 (cl)) / a2 (D3) 

i - cos' 1 (cos (cl) / c2) (D4) 

where i is the orbit inclination in degrees. Equations (Dl) and (D2) are 
iterative steps used to calculate the orbit inclination. 

The following equations used to calculate orbital parameters give the apparent 
coordinates of the Sun to a precision of 0.01 degrees and time to a precision 
of 0.1 minutes. Similar equations can be found in references 3 and 6. As 
expected, the location of the Sun is dependent on the day of the year. Thus a 
time argument is used to calculate the Sun's location in a right ascension - 
declination coordinate system. The time argument is expressed as (ref. 5) 

n =- -4018 . 5 + d (D5) 

where d is the day of the year. The time argument is used to calculate 
several orbital angles 


L - 280,460 + 0.9856474 n (D6) 

g * 357.528 + 0.9856003 n (D7) 

X “ L + 1,915 sin g + 0.02 sin 2g (D8) 

e - 23.439 * 0.0000004 n (D9) 


where L is the mean longitude of the Sun, g is the mean anomaly, X is the 
ecliptic longitude, and €. is the obliquity of the Earth's equator to the 
ecliptic, in degrees. These equations are used to calculate the right 
ascension of the Sun (D10) , and the solar declination (Dll) . 


u 

- tan*^ (cos € 

tan X ) 

(D10) 

§ 

- sin‘^ (sin €■ 

sin X) 

(Dll) 


In equations (D6) through (Dll), and other equations involving angular results 
in degrees, the results often produce angles not within the customary range of 
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0 to 360 degrees. Therefore, all angles should be converted to values within 
this range. Using the results from equations (D8) and (D9) the Sun's 
position, normalized with respect to the distance from the Sun to the Earth, 
can be expressed as a vector with rectangular coordinates. 


X g - cos X 

( D 1 2 ) 

Y g - cos £ sin X 

(013) 

Z s — sin sin X 

(014) 


The solar vector in the right ascension - declination coordinate system is 

S - X s i + Y s j + Z s k (D15 ) 

The spacecraft position can be expressed in a similar manner using equations 
dependent on the day of the year, sun- synchronous orbit time, and the 
inclination angle. The time argument for these equations is (ref. 4) 

T - ( J - 2415020 ) / 36525 (D16) 

where T is the time in Julian centuries from the date January 0.5, 1900 and J 
is the Julian date. The right ascension of the Greenwich meridian at midnight 
can be expressed as 

r - 99.6909833 + 36000.7689 T + 0.00038708 T 2 (D17) 

with the result falling between 0 and 360 degrees. The equation for longitude 
of ascending node is 

- r + 15.0410682 t 

where t is the sun- synchronous orbit time in hours past midnight. 

(D4) and (D18) can be used to calculate the spacecraft position in 
coordinates for a right ascension - declination coordinate system, 
coordinates are 


(D18) 

Equations 

rectangular 

These 


- sin i sin 

(019) 

Y n - - sin i cos -ru 

(D20) 

Z n = cos i 

(021) 

and the corresponding vector equation 

N “ i + Y n j + Z n k 

(D22) 


The angle between the orbit plane of the spacecraft and the solar vector can 
be found by taking the dot product of the two vectors in equations (D15) and 
(D22) . This result yields the "Beta Angle", in degrees, which is one of the 
most important parameters for spacecraft thermal analysis. 

p - 90 - cos’ 1 ( S • N ) (D23) 
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